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A constant in the expression t = (see Ref. 5) 

f[(Y-l)/2) H//(yTg))’'“ 

C((y- 1)/21 MeV(yT^))''^ 

t(l+[(Y-l)/2] y)/(T^/Tg)}-l 
((l+(ft)*'“ [(y- 1)/2] Ke2)/(T„/Tg))-1 

a constant in Lav/ of the Wall (usually 5.1) 

Skin friction coefficient ■ry(V2)pg 

5.1 - 0.614/aK + (1/K) In («u^/v^.p 

Constant in mixing length (usually 0.4) 

Mach number 
Pressure 

Turbulent Prandtl number 
Temperature 

Velocity in streamwise direction 

(Ug/A) arc sin {[(2A^u/Ug) - B]/(B^+4A^)'^> 

(Ug/A^) arc sin {[{2h^^u/u^) - B^ 1 /(B^ 2 + 

Friction velocity (t„^/p^^)'^ 

u* + sin (3/(3^'’t-4A^)'^} 

Coles' universal wake function 
Coordinate noraal to v/all 
Ratio of specific heat 
Boundary layer thickness 
Displacement thickness 
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Momentum thickness 

y/6 

Kinematic viscosity 
Coefficient of Wake Function 
Density 
Shear Stress 


Subscripts 


e 

0 


V/ 


Freestream conditions 
Stagnation conditions 
Conditions at the wall 



INTRODUCTION 


The wall-v/ake velocity profile has been used to represent turbulent boundsr 
layer profiles for both adiabatic and non-ad1abat1c flows and for flov/s with 
or without pressure gradients [1-5]. A least squares fit of a wall -wake velo- 
city profile to an experimental velocity profile may be used to determine 
and 6 for the profile. An accurate representation of the mean velocity distri- 
bution in the turbulent boundary layer can be very useful in integral analyses 
of turbulent flow problems; in the analysis of flows in which strong interactions 
occur (as for example in shock wave-boundary layer interactions) and in which com- 
bined viscid-inviscid analyses are required, a profile which provides a good re- 
presentation of both and 6 can be quite important. With most earlier versions 
of the wall -wake profile the velocity gradient at the boundary layer edge has a 
non-zero value. The values of 6 for these profiles may, as a consequence, cor- 
respond to points in the flow where viscous effects are substantial, especially 
in very high-speed flov/s. The purpose of this note is to suggest a form of the 
wall-wake profile which is applicable to flows with heat transfer, and for which 
9u/ay - 0 at y =6. The modified profile, which takes into account the effect of 
turbulent Prandtl number, has been found to provide a good representation of 
experimental data for a wide range of Mach numbers and heat transfer. The 
values v/hich are determined by a least squares fit of the profile to the data 
agree v/ell with values which v/ere measured by the floating element technique. 

In addition, the values of 6 determined by the fit correspond more closely to the 
outer edge of the viscous flow region than those obtained with earlier versions of 
the wall-v/ake profile. 
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VELOCITY PROFILE 


The difference between the profile to be discussed here and some of the 
earlier versions is due to differences in the lav/-of-the-wall component. 

Van Driest [6] developed a compressible law-of-the-v/all under the assumption of 
a turbulent Prandtl number of unity and a constant shear stress near the wall. 
It may be v/ritten as, 
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With the addition of a wake component, Eq. (1) becomes. 
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The left-hand side of this expression is the transformed velocity, U , v/hich 
is used in References 3 and 4. 

The use of Eq. (1) in a v/all-wake profile results in a non-zero value of 
the velocity gradient at the boundary layer edge. If we assume, as in Reference 
5, that T may be expressed as t = 'r^^{l"n^)> and if for Pr^ f 1 the temperature 

distributions may be v/ritten as (T/T^P = 1 + "* ^ (cf. Schlichting 

[7]), with (p/p^) = 1/(T/T^), we may follow the procedure used in Reference 5 
to obtain a v/all-wake boundary layer velocity profile of the following form. 
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For a ®, v/hich corresponds to the assumption of a constant shear stress dis- 
tribution in the derivation of the compressible law-of-the-wall , and for 
Pr^ -j- 1, u** u*, as in Eq. (2). Setting n = 1 in Eq. (3) yields the follow- 


ing expression for n/K, 
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For mathematical convenience [1-3] we may replace W(n) in Eq. (3) by 
cos hit] and v/rite the wall-v/ake velocity profile as. 
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The method of least squares has been used to fit the wall -wake profile, 

Eq. (5), to experimental velocity profiles reported by Hopkins and Keener [8], 
Voisinet and Lee [9], Horstman and Owen [10] and Kilburg [11] for zero pressure 
gradient flov/s with heat transfer. The walV-v/ake velocities were obtained under 
the assumption of Pr^ = 0.8 and a = 1, the latter corresponding to the assump- 
tion of a linear shear stress distribution in the derivation of the law-of-the- 
wall. K was taken as 0.4 and C and 5.1. Data in the sublayer were excluded. 



4 


The results shovm here are from studies for which directly measured values of 
wall shear stress were reported [8, 9, 10]. The data reported by Voisinet and 
Lee and Hopkins and Keener were obtained on wind tunnel walls. The Horstman- 
Ov/en data are for flow along the cylindrical surface of a cone-cylinder model 
whose axis was aligned with the primary wind tunnel flow direction. 

Comparisons of the experimental and wall -wake profiles are given in 
Figure 1. Also shovm on the figures are values of 6*, 6, C|r and ll/K as deter- 
mined with the wall -v/ake profile. Experimental values of 6*. 0 and as ob- 
tained from the references are also listed. The wall -wake profile is seen to 
provide a good representation of the experimental velocity profiles. In addi- 
tion, the values determined by the fit of the proposed profile agree rea- 
sonably well with the directly measured values. For the data of References 8 
and 10, the values of 5 determined by the fit correspond to points in the flow 

at which u/u , Pq/Pq and T^/T^ are essentially unity. If a is set equal to 

0 0 

« in Eq. (5) the values of 6 are on the order of five percent lower than for 

a = 1 and the values of u/u^, Pq/P^ and 6 are correspondingly lower. For rea- 

e 

sons which are not known at this time the value of 5 determined for the data of 

Reference 9 corresponds to a point in the flow at which u/u^ = 0.965, Pq/Pq = 

0.56 and T^/Tq = 0.938, even though the wall -wake profile and the data agree 
e 

very closely between the wall and the value of 6 determined by the least squares 
fit. The authors comment on the temperature and history effects in some of 
their profiles. The rather Tow value of 6 may be associated with these effects. 
For a -x^t 6 is about five percent lov/er than for a = 1. 

The wall-wake values of 5* and 8 agree well with expertmental values for 
the data of Reference 10. The larger differences for 5* and 9 for the data of 
Reference 8 are due to differences between the temperature distribution for the 
experimental profile and the distribution for the wall-v/ake profile with 
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Pr^ - 0.8. In the case of the profile from Reference 9, the larger differences 
appear to result from the fact that the velocity, total temperature and total 
pressure continue to increase v/ell beyono the value of «s determined by the 
least squares fit. 

The values of n/K shown in Figure 1 were determined from Eq. (4), assum- 
ing Pr^ = 0.8 and a = 1 and v/ith u^ and 6 determined from the wall -wake fit. 

Note that the expression for n/K differs from that used in References 2-4 
through the inclusion of the Prandtl number effect, through the term 0.614/aK, 
and through the differences in u and 6 which are obtained with the modified 

T 

profile. Although the pressure gradient was zero for the three profiles shown, 
the values of n/K are quite different. The values varied considerably from profile to 
profile for all of the data examined and tended generally to be higher than for 
low-speed or compressible adiabatic flows. 

The experimental profiles may also be plotted in the form (u** - u**)/u^ 
versus y/A, where A is analogous to the defect thickness used by Gran et al., 

[ 4 ] in their comparison of cold-wall, high-speed data with low-speed data [13]. 

Such a plot is shown in Figure 2, While the profile from Reference 9 deviates 
from the low-speed results, the profiles from References 8 and 10 correspond 
quite closely with the low-speed data. This suggests that the latter tv/o pro- 
files were nearer equilibrium than the former, v/hich, in turn, may account for 
the fact that the value of 6 determined for the profile from Reference 9 does 
not correspond as closely to the viscid-inviscid flow boundary as was found for 
the other two profiles. 
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Figure 1 
Figure 2 


Velocity distributions in boundary layer v;ith heat transfer. 

Cornparison of generalized velocity defect for lov/-speed flow 
and compressible flow v/ith heat transfer. 
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